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We present an ab initio study of the ground state of an ideal coupled two-component gas of 
ultracold atoms in a one dimensional optical lattice, either bosons or fermions. Due to the internal 
two-level structure of the atoms, the Brillouin zone is twice as large as imposed by the periodicity 
of the lattice potential. This is reflected in the Bloch dispersion curves, where the energy bands 
regularly possess several local minima. As a consequence, when the system parameters are tuned 
across a resonance condition, a non-zero temperature phase transition occurs which arises from an 
interplay between internal and kinetic atomic energies. For fermions, this phase transition is of 
topological character since the structure of the Fermi surface is changed across the critical value. 
It is shown that these phenomena are also expected to occur for two and three dimensional optical 
lattices. 

PACS numbers: 37.10. Jk,05.30. Jp,05.30.Fk,03.75.Mn 
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I. INTRODUCTION 



The coupled dynamics between an ultracold atomic 
gas and an optical lattice has gained enormous attention 
during the last decade. Originally, D. Jaksch et al. pro- 
posed that the Mott- super fluid quantum phase transition 
(PT) can be realized using cold atoms dipole coupled to 
a standing wave laser field A mere four years later, 
this PT was observed in the seminal experiment by I. 
Bloch and coworkers One among the many reasons 
for the great interest in these systems is the possibility 
to experimentally realize, and therefore verify, theoret- 
ical models developed in the field of condensed matter 
physics, see the review Several achievements in this 
field have been accomplished both theoretically and ex- 
perimentally, and to only mention a few; disordered sys- 
tems and matter localization [3], frustrated systems Q, 
the strong atom-atom interaction regime and the Tonks 
gas [sl, cold atom quantum Hall counterp arts 7] and 
quantum information processing [8] . 

Almost exclusively, effective dispersive models have 
been considered, where the internal level structure of 
the atoms can be discarded and one is left with a single 
atomic level. Exceptions are works on spin-dependent 
atoms, which have shown to acquire new phases com- 
pared to the spinless case [9]. Another, related system 
is mixtures of different atomic species, again containing 
novel physics However, these models do not analyze 
direct coupling between the internal atomic states, as will 
be the topic of the present paper. Furthermore, atom- 
atom scattering plays a crucial role in the dynamics of all 
above references. For example, the Mott-superfluid PT 
in the Bose-Hubbard model arises from a competition be- 
tween the atomic kinetic energy (hopping between neigh- 
boring sites) and the onsite scattering interaction [ll|. 
Here we will show that by taking into account for dipole 
induced transitions between coupled atomic states, PTs 
can occur even f or ideal gases lacking atom-atom inter- 
action. This originates from the competition between the 
atomic kinetic and the internal energies, in comparison 



with kinetic and scattering energies. 

Consequently, we do not restrict the character of the 
atom-field interaction to be dispersive. In fact, we al- 
low for a vanishing atom-field detuning and in particular 
study the behavior for both positive and negative detun- 
ings. Coupling between internal levels of atoms in optical 
lattices is rather unexplored. The first analysis seem to 
date back to the works by Krutitsky et. al., where they 
considered effective coupled two or three level models de- 
rived from Raman interactions [l^, [l^ . They found that 
the hopping term in the corresponding Bose-Hubbard 
model may be tuned to be either positive or negative, 
leading to new phases. An internal coupled model was 
as well studied by Garcia- RipoU et. al. in Ref. ^1^. 
Contrary to [H, [3 , we consider here a direct coupling 
between two atomic levels, which is actually what one 
gets in the standard optical lattice model by decreasing 
the amplitude of the detuning. Another crucial differ- 
ence between the current contribution and the ones of 
Ref. [H, is that our model goes beyond any tight 
binding approximation, in which hoppinig beyond near- 
est neighbours is allowed. The detuning in our model 
may be seen as having similar role as the parameter 9 
(relative phase between the two Raman lasers) in [l^ . 
An important point to note is that in Ref. T4], on the 
other hand, the coupling between internal levels is not 
originating from a non-zero detuning (as in the current 
paper) or from a non-zero angle 9 (as in [T^ . [l3|). but 
from atom-atom scattering. We explicitly show that the 
two-level structure of our model renders dispersion curves 
having local minima. It is known that such phenomena 
can give rise to topological PTs [H, [H, which in- 
deed are found in our model as well. A topological PT 
(for fermions; a topological change in the Fermi surface) 
does take place when the detuning changes sign. The 
same kind of PT indeed occurs also for bosons. The 
non-zero temperature situation is also considered and the 
topological PT is found to be stable against temperature 
fluctuations. 

The outline of this paper is as follows. We first discuss 
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the single particle Hamiltonian in Sec. [TTl and point out 
some of the symmetries associated with it. One of these 
symmetry operators indicates that the Brillouin zone ex- 
tends over twice the size expected from the periodicity of 
the optical lattice, a fact that is clarified even further in 
the proceeding Sec. IIIII In Sec. IIIII we begin by analyz- 
ing the energy spectrum and the two lowest band Wan- 
nier functions. Using this knowledge we demonstrate the 
presence of a PT, both for fermions and bosons at zero 
and non-zero temperatures. We also discuss the effect 
of weak atom-atom interaction. A discussion of possible 
extensions is left for the conclusions given in IIVI Last, 
in the appendix IIVI we also present an effective Raman 
coupled model that would provide qualitatively the same 
results, but with the benefit of using two meta-stable 
atomic states. 



II. SINGLE PARTICLE HAMILTONIAN 

The Hamiltonian for a single two-level atom, whose 
internal states dipole couple via a standing wave laser 
field, reads 

jp hA ^ 

H ^ \- —az + 2hgcos{kx)aoo- (1) 

2m 2 

Here, p and x are atomic center-of-mass momentum and 
position, m its mass, A the atom-field detuning, g the 
effective atom-field coupling and k the field wave num- 
ber. The internal states of the atom are labeled |±) 
and the Pauli matrices operate as ^zji) = and 
(T2;|±) = 1^). Before proceeding we introduce dimension- 
less variables through the characteristic length and 
energy Er = h'^k'^/2m; 



As a periodic operator, the eigenstates and eigenvalues 
of Hamiltonian ([3]) are characterized by two quantum 
numbers, one discrete band index v and a continuous 
quasi momentum q 



v=l,2, 3, 



(5) 



The range of q assigns the first Brillouin zone and is de- 
termined by the symmetries of the problem. For any 
parameters, the operator T — e^^^ commutes with the 
Hamiltonian, where A = 27r is the field wavelength in 
scaled units. Not as evident is that also the operator 



(6) 



is a symmetry of the Hamiltonian [20|, |2l|. Thus, there 
is a natural tt periodic structure of the Hamiltonian. 
This signals that the first Brillouin zone extends beyond 
— l/2<g<l/2as implied by the 2-k periodicity of the 
T operator, and implicitly we have 



(-1,1] 



(7) 



However, for A = the spectrum is doubly degenerate 
and in this limit it might be more convenient to define 
the Brillouin zone within — 1/2 < g < 1/2. One may 
note that T = . The physical background of this addi- 
tional symmetry operator / emerges from the fact that 
absorption or emission of a single photon flips the internal 
states, |±) |=f), while only for a two-photon process is 
the internal state unchanged. 



III. STRUCTURE OF THE GROUND STATE 



Model characteristics 



= ka 



A - 

Er 



Tig 



(2) 



In the 



and 



becomes in scaled variables 



nomenclature, Eq. ([T]) 



H = - 



A 
"2" 



2g cos(i) 
A 



2gcos{x) - — 



(3) 



which serves as our model Hamiltonian. We note that 
the external field couples the hare states |±) and simul- 
taneously shift the momentum by either ±1. For A = 0, 
the unitary operator U — ^ [a^ + CTz) decouples the in- 
ternal levels, and one obtains two Mathieu equations [l^ 
with diahatic potentials V^{x) — ±2gcos{x) [l9j . The 
internal atomic states of the decoupled equations are 



|l> = ?^l+> = ^(l+> + |-», 



(4) 



|2)=t/h 



T2' 



To find the ground state of N non-interacting two- 
level atoms we need to diagonalize ([3]); both within the 
internal and the motional degrees of freedom. We first 
note that for a given quasi momentum q, the Hamiltonian 
can be written on the block form H — (E) H^, and the 
states coupled by the sub-Hamiltonians read 



-v)\- 
-v)\- 



\q + v)\+) 
\q + v)\-) 



rj even 
77 odd 

rj even 
77 odd. 



(8) 



where the first ket is the momentum eigenstate p\q + ri) — 
{q + Ti)\q + rj) and 77 is an integer corresponding to the 
momentum shift governed by absorption and emittance 
of photons. This decoupling of the dynamics is of great 
importance as a given number of atoms residing in each 
subset is constant when the parameters are varied. The 
states ([5]) are eigenstates of the Hamiltonian in the ab- 
sence of matter-field coupling 5 = 0, with bare energies 



(9) 
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where the ±-sign is different for the two sets ([8]). In 
this trivial situation, the detuning A simply shifts the 
parabolic dispersions. The general eigenstate, Block 
state, were introduced in Eq. ([5]). The eigenvalue E'^{q) 
is the i^'th energy band/dispersion curve, given for q e 
(—1,1]. The eigenstate in a position representation can 
be written using its constituent Block states as 



(10) 



Note that the above specifies just two possible con- 
stituent Bloch states and any other internal basis would 
define two new ones. In particular, ipl''^{x) and ip^y{x) 

are related via the operator [/ = {(Jx + (Tz) of Eq. Q. 
In terms of the the dispersions, the Hamiltonian is given 
by 



H 



E E 

2^=1 ge(-i4] 



(11) 



where h'^ is the operator giving the number of atoms with 
energy (q) and here we have introduced a chemical po- 
tential /i. Note that for a given atom number N one has 



95 (■ 



N, where (n^) is the expectation 
value for the particular state of interest. The Bloch states 
can be expressed in the bare basis (|5]) 



\M<i)) = { 



E <(9)l^.('z)> 

7— — oo 

oo 

E d';{q)\M<l)) 



(12) 



k »?=-oo 

Note that expanding \ip^{q)) does not mix the states 
\ip^{q)) and |</',,(g)) due to the block-diagonal form of 
the Hamiltonian H = H^p O i?^ in this bare basis. We 
emphasized that a quasi momentum eigenstate, or Bloch 
state, is a linear combination of the two internal atomic 
states |±). The same holds for the Wannier functions 



= w+{x-R)\+)+w~{x-R)\-) 
= wlix^R)\l)+wUx-R)\2). 



(13) 



As defined above, the constituent Wannier func- 
tions w'^ {x — R) and w'J (x — R) are normalized as 
J2i=± I I'^li^ ~ R)\^dx = 1, and likewise for wl{x — R) 
and w^x — R). Like for the Bloch states, the constituent 
parts depend on the internal basis. Equipped with this 
machinery, we now turn to analyze the complex band 
structure of the Hamiltonian ([3]) . 

Figure [T] displays some examples of the band structure 
- the dispersion curves, where the first three bands are 
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FIG. 1: The three lowest energy bands of Hamiltonian ([Sj, 
where the scaled dimensionless parameters are g — 0.1 and 
A = 0, 0.25, 0.5, 1 in plots (a)-(d) respectively. The two low- 
est bands are distinguished between the two sub-Hamiltonians 
Hip (solid line) and (dashed line). 



shown. The resonance case, A = 0, is given in (a) and 
here the energy spectrum is doubly degenerate as pointed 
out in the previous section. However, for non-zero A, 
the degeneracy is lifted (indicated in (b)-(d)) and as a 
consequence, the first Brillouin zone extends over quasi 
momenta q € (—1,1]. The dashed line corresponds to 
the lowest dispersion curve of the sub-Hamiltonian H^, 
while the other lowest line derives from H^p instead. Note 
that the dispersion curves have more than a single mini- 
mum in all plots. The non-monotonicity of the dispersion 
curves can be understood from the interplay between the 
different energies; kinetic energy causing the parabolic 
structures, the coupling energy which typically splits the 
degeneracy and the internal atomic energy that (loosely 
speaking) shifts the dispersions. 

Let us briefly discuss the dispersion curves and their 
limiting cases in more detail. In the regular dispersive 
situation one has A ^ g and the coupling between the 
internal adiabatic states can be neglected. To verify this, 
let us introduce the adiabatic states p3 ] as columns of 
the unitary operator 



Uad 



cos(6l/2) sin(6l/2) 
sin(6'/2) -cos(6'/2) 



where 



tan(6') 



lcos{x) 
~A ■ 



(14) 



(15) 



The operator Uad then diagonalizes the 2x2 matrix of 
the Hamiltonian ([3]). However, due to its a;-dependence 
will it not commute with the momentum operator p. This 
causes non-diagonal terms in the transformed Hamilto- 
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nian [22 



(16) 



Here, Vcent(x) is a centrifugal term turning up on the 
diagonal [Ij], ri(a;,p) is the non-adiabatic coupling 
and specifically 



^±"'(^) = ±y (y) +Vcos2(x) (17) 

are the adiabatic potentials. It follows that the centrifugal 
and adiabatic correction terms are small in the A ^ g 
regime [13], giving the adiabatic potentials 



yad 



^ A _L VC0S2(X) 
(x)-±y± ^ . 



(18) 



Thus, we derive the regular situation most commonly 
considered in the literature and we especially note that 
the first Brillouin zone extends over q E (—1, 1]. Choos- 
ing A > and considering the weak coupling limit 
Ag'^/A 0, the states I?)]— ) with q e (—1, 1] have eigen- 
values 



E'iq)=e'=q'- 



A 



(19) 



For the states \q ± 1) 
dispersions 



on the other hand, we have the 



i?"(g)-(9±i) 



A 
2"' 



(20) 



However, as is noticeable from Eq. the states \q)\—) 
and |g± 1)|— ) are not coupled and when one restricts the 
analysis to just one of these sets and one consequently 
regains a regular spectrum, by which we mean disper- 
sion curves that posses only a single minimum within 
one Brillouin length, dE^{q)/dq > for < g < 1. More 
interesting is the intermediate regime, neither adiabatic 
nor diabatic (A = 0), where the above states may indeed 
coexist if we assume that both of the states \q)\±) can 
be present. Thus, when the couphng between |±) states 
cannot be neglected, one must take into account all the 
states and the lowest energy band contains minima at 
q = ±1 as well as for q — 0. 

As pointed out, here we allow for atomic states |-|-) 
and |— ) to have the same momentum. It is understood 
that the same arguments hold for A < 0, making the 
replacement |+) ^ |— ). We should mention that even if 
we only consider atoms within one subset ([5]) , the disper- 
sion curves will contain several local minima 21, 24], and 
the results presented are valid also in such cases. In the 
diabatic limit, A = 0, we saw that we can separate the 
dynamics into two uncoupled problems with diabatic po- 
tentials V^{x) — zL2gcos{x). The spectrum is then dou- 
bly degenerate and the first Brillouin zone is therefore 



most properly defined within q G (—1/2, 1/2]. Notewor- 
thy is the fact that the diabatic potentials have minima 
either for x'},^ = 2mT or = (2n + 1)tt for integer n, 
while the minima for the adiabatic potentials are either 
x'^ — nn or xj^ = + l) ±-sign we indicate 

the corresponding diabatic or adiabatic potential. 



R=0 



R=7i/2 




FIG. 2: The constituent squared amplitudes of the Wannier 
functions of the first, (a) and (b), and second, (c) and (d), 
Bloch band. The left plots display the Wannier functions for 
R = 0, while the in the right plots R = n/2, corresponding to 
the minima of V!^'''''^{x) and V+'*(a;) respectively. Solid lines 
depict \wf{x — _R)p and dashed line (x — R)\^- The other 
dimensionless parameters are g = 0.1 and A = 0.25. 



The coupled two- level character of the system gives rise 
to rather peculiar constituent Wannier functions {x — 
R) or wl''^{x — R). In general, R is chosen such that the 
potential is minimal at a; = i?. However, in the present 
model R is not a priori given since the potential has a 
complex two-level structure. From the symmetry of the 
Hamiltonian it follows that {x — R) ^ w~ {x — R) 
for A <-> —A, while \wl,''^{x — R)\ are invariant under 
such sign change. In Fig. [5] (a)-(d) we show examples of 
\wf{x - i?)]2 (a) and (b) and \wf{x - R)\^ (c) and (d), 
for i? = or i? = 7r/2. The dashed line shows \wY {x — 
i?)]2 and the solid \w^{x - R)\'^. Note that for R = 
0, the two constituent Wannier functions are identical, 
which, however, does not hold for R = Tr/2. For R = 0, 
the constituent Wannier functions resemble typical ones 
obtained in one component systems [25j . This is not tru 
e for R — Tr/2 (corresponding to the minima of V_^'''{x)), 
where only \w~{x — i?)]^ shows the regular shape and 
{x — i?)]^ looks as if R coincide with a maximum of 
the potential. Note that the result of the figure presents 
an intermediate regime where neither the adiabatic nor 
the diabatic approximations can be considered; g = 0.1 
and A = 0.3. 
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FIG. 3: Coefficients c^,((?) (a) and d^(g) (b) defined in 
Eq. (|12|l . The parameters are as in Fig. [T] (b), but with a 
larger coupfing g = 0.5. 



B. Zero temperature many-body ground state 

We will focus on the lowest band, v — \, thus assuming 
a filling factor ^ = ^ < 1 , where N is the total number of 
atoms and K number of sites. In the numerics, K will be 
taken large enough (typically K > 100) to assure small 
boundary effects. The N many-body ground states for 
fermions and bosons, at zero temperature and for filling 
factors ^ < 1, are 

i*>F=n4io)^nK): (21) 

geS qeQ 

|*)b = |no,n+i), no + n+i=N (22) 

respectively and Q contains those values of q which are 
inside the Fermi sea, is the Fermi creation operator 
of mode q of the lowest Bloch band, |0) is the vacuum 
and Hq (= (hq)) again characterizes the number of atoms 
in quasi momentum mode q; X^ggg"-? ~ Note that 



for bosons, the ground state is degenerate and no and 
n+i may pertain any positive values such that the total 
number gives (we have not included q — ~1 as it lies 
outside the Brillouin zone), and in fact any linear com- 
bination of these degenerate states is adequate. Further, 
the state ([22]) is valid for any filling factors In fact, 
the same holds in general for any non-integer filling fac- 
tor ^. As states (PT|) and are prodcut states, there 
exist no correlation between the atoms. Nonetheless, en- 
tanglement between atomic motion and internal atomic 
states occurs for each atom (provided A, g 0). This is 
clear since 

{x\Mq))^x{^)ia\+)+b\-)), (23) 

for some normalized function x(a;) and |ap + |5p — 1. 
The expansion of the Bloch eigenstates in terms of bare 
states is given in Eq. (fT^ . The corresponding coefficients 
are displayed in Fig. [3l for the lowest Bloch band with 
parameters g — 0.5 and A ~ 0.25. For \q\ < 1/2, the 
c^(g) are all zero, while djj^{q) = for 1/2 < |g| < 1. 
A result deriving from the block structure of the Hamil- 
tonian ([8]). The cj^{q) coefficients are dominated by the 
ones with either ^=lor/i = — 1, and = is the 
most prominent coefficient of the dj^{q). Note that the 
probability for the atom to be found in the state |— ) for 
a given eigenstate ip,^{q) is 

Pi-;Mq))= E E l<(9)P (24) 

V even jj odd 

and P{+;Mq)) = l-Pi-,MQ))- 

We define the internal density per particle as 

'^^(^) = ^EEl^^.(^)l'' (25) 

where the second sum runs over occupied momentum 
states, and N is the number of atoms. The density p^{x) 
for fermions is shown in Fig. [?] (a) as function of the 
detuning A where the lowest Bloch band is half filled, 
^ — 1/2. For other filling factors one regains very similar 
plots. By substituting A ^ — A, the p^{x) is identical to 
p~{x). Noticeable is the discontinuity at A = 0. Thus, 
around resonance, A = 0, the population of |-|-) and |— ) 
atoms may fluctuate and in particular a first order FT in 
the collective atomic inversion (per particle) 

^ ^E E [P{+;Mq))-PHM<i))] (26) 

is expected. Note that the above inversion ((26)) may be 
derived from the expectation value of the collective in- 
version 



6 




-2-1 12 

A 

FIG. 4: The upper plot (a) displays the Fermi particle density 
p~{x) (|25p for the internal state | — ) at half filling of the lowest 
Bloch band. The lower plot (b) shows three examples of the 
collective atomic inversion jz (|26p as function of A for bosons 
(dashed lines) and for fermions (solid lines). In (a) g — 1. 



where az is the i'th atoms z PauU matrix measuring 
the i'th particle inversion. In the Dicke model of N two- 
level atoms interacting with a quantized field mode, jz 
is often considered as an order parameter having a dis- 
continuous first order derivative at the critical atom-field 
coupling which defines the normal- superradient quantum 
PT [26|. The inversion jz as function of A, presented 
in Fig. 13] (b), illustrates the discontinuity incorporated 
in our model. Note that this PT is different from the 
one of the Dicke model, where the PT originates from 
competition between field and interaction energies and 
not between kinetic and internal atomic energies as in 
this paper. In fact, the Dicke model lacks any sort of 
periodic lattice potential. 

We have verified that the PT obtained in our model 
also occurs when one restricts the dynamics to one of the 
subsets Assuming a single particle in either of the 
subsets dH), one has -Bi(±l) = Ei{0) for A = 0, while 
i?i(±l) 7^ Ei{0) for A 7^ (clear for example from the 
dashed curve if Fig. [1]) . Note that this is not true if we 
include both subsets where i?i(±l) = £'i(0) is always 
true (here Ei (q) is the lowest Bloch band of either iJ^ or 



H^). Thus, the single particle ground state energy is ei- 
ther Ei{±l) or Ei{0) depending on the sign of A, and as 
one tunes A across resonance the atom will absorb/emit 
one photon and W ^ — W. 

The same argument holds also for many atoms where 
each one of them swaps internal states, |±) ^ |=f), while 
passing through A = 0. For a fermionic system, the 
PT is of topological character, since the structure of the 
Fermi surface changes across the critical point due to the 
fact that each atom emits/absorbs a momentum kick. 
The mechanism behind the PT is the interplay between 
internal and external atomic energies. The collective in- 
version, related to internal structure, reveals the PT, but 
its topological character is, however, only manifested in 
the external degrees of freedom. Even for bosons can 
the PT be termed topological becuase the nature of the 
atomic phase space distribtion is altered while crossing 
the critical point. Topological PTs, or Lifshitz transition, 
have been studied comprehensively in quantum Hall and 
superconducting systems [H, [l^. However, a recent pa- 
per considered a topological PT of ultracold fermionic 
atoms in an anisotropic three dimensional optical lattice 

The non-linearity, imposed by including atom-atom in- 
teraction, is known to cause changes such as loops in the 
band structure around the avoided crossings of the dis- 
persions [13] ■ This effect occurs for rather strong cou- 
plings, or more precisely when the effective scattering 
amplitude becomes comparable or larger than the lattice 
depth [27|. For reasonable temperatures is the Fermi 
surface well below the first band gap where the greatest 
effects of the non-linearity is expected to manifest itself. 
Consequently, the topological PT is supposed to be in- 
sensitive for such modifications, since the general struc- 
ture of the spectrum remains and thus the topological 
change in the Fermi sea would still take place as the de- 
tuning is tuned across resonance. Nonetheless, one may 
expect that even for weak interactions, the interaction 
induced splitting of the degeneracy of the ground state 
might cause a re-formation of the ground state into a 
correlated state which is fundamentally different from the 
simple product state. Thus, the ground state is no longer 
a simple product state of Fock states as given in Eqs. (pij) 
and (|22p. On the other hand, the inversion is invari- 
ant for such re-configuration of the ground state because 

{Mq)\'^z\Mq)) = (V^.(-g)|a,|v^,(-q)), (-1,1]. 

One limitation of our model is spontaneous emission 
of the excited atomic state. This can be circumvented 
by considering instead a Raman coupled A-atom setup, 
in which the excited state is adiabatically eliminated due 
to dispersive interactions. The resulting effective Hamil- 
tonian, derived in the appendix, has a similar form to 
the one of Eq. ([3]) apart from two additional Stark shift 
terms. By proper choices of parameters can these two 
terms be made small and the corresponding spectrum 
then possesses the same properties as Hamiltonian ([3]). 
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C. Non-zero temperature ground state 
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0.1 0.2 0.3 0.4 0.5 

T 

FIG. 5: The collective atomic inversion jz (a) for fermions 
(black lines) and bosons (gray lines) as function of the de- 
tuning A and different temperatures, fermions: T = (solid 
line), T = 0.02 (dotted) and T = 0.1 (dot-dashed) and bosons: 
T — (solid line) and T = 0.1 (dashed line). The inset is a 
close-up of the curves to the right of A = 0. The second fig- 
ure (b) depicts the temperature dependence of the inversion 
gap Sjz of Eq. 1281 for Fermions (solid black line) and bosons 
(dashed gray line). In both plots, g = 1/2 and the chemical 
potential is chosen such that the number of particles is half 
the lattice number. 

We turn now to the situation of non-zero temperature 
and thermal excitations. As argued above, each eigen- 
state \'ipi{q)) in the zero temperature ground state un- 
dergoes a coUective Rabi flip by sweeping the detuning 
across resonance, which suggests that also for excitations 
around the Fermi surface or of the Boson ground state 
will the atomic inversion be discontinuous around A = 0. 
Let us introduce the gap function 

Sjz = lim - hm (28) 

and study the behavior of 6jz for various temperatures T. 
The temperature dependencies of the collective atomic 
inversion jz and the inversion gap Sj^ are presented in 
Figs. [H] (a) and (b) respectively. The gap is slightly larger 
for Bosons than for Fermions. Surprisingly, the gap 5jz 
approaches a non-zero value (~ 0.75 for this choice of 



coupling; g = 1/2) in the large temperature limit. Even 
at fairly large couplings, this asymptotic value is non- 
zero. 
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FIG. 6: Population of excited bands Pex, in the case of 
fermions, as function of the detuning and of temperature; 
T = 0.02 (solid line), T = 0.1 (dashed line) and T = 0.2 
(dot-dashed line). At resonance and moderate temperatures, 
the excited population approximates 1/2 due to the degener- 
acy at that point. For large detuning, the gap between the 
two lowest bands increases causing a decrease in the excited 
band population. The inset shows a close-up of the popu- 
lation around A = 0, indicating that for large temperatures 
Pex exceeds 1/2 at resonance due to occupation of third and 
higher bands. Here g = 1/2 and ^ — 1/2. 



Another relevant question is the amount of populations 
of excited quasi momentum states \'ip^{q)) for non-zero 
temperatures, especially for fermions. The two lowest 
Bloch bands are identical at resonance, A = 0, and conse- 
quently equally populated. For small temperatures, only 
states corresponding to these two bands are substantially 
occupied, while for large temperatures also the third and 
forth band will be populated. For non-zero detuning, the 
degeneracy is lifted and the populations of the two low- 
est bands are no longer balanced. We may define the 
individual band population 

E K)' (29) 

where (np = Tr [pT?^g] is the expectation value of the 
number operator h'^ for the state pr given at temperature 
T. Thus, Pi, measures the population in band v, and 
in particular for A = do we have Pi = Pi+i for i = 
1, 3, 5, ... due to the degeneracy. The total population 
of excited bands is given by 

oo 

Note that P^x includes the second band v — 2 which for 
A = is degenerate with the lowest band v = 1 and in 
this special case is the subscript ex (excited) misleading. 
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Naturally, < Pex < I and also 1/2 < Pex(A = 0) < 1. 
Figure [5] displays Pex as function of A and for three dif- 
ferent temperatures. As expected, the excitations in- 
crease for a large temperature and decrease for a large 
detuning. The inset gives the population close to reso- 
nance, and it is in particular seen that Pea; (A = 0) > 1/2 
(Pea; (A = 0) = 0.5025 for T = 0.2) from the fact that 
the bands = 3, 4, ... begin to be populated. The shape 
of Pex seems fairly Lorentzian, but the inset reveals that 
Pea; is indeed not Lorentzian in the vicinity of A — 0. 
Again we restrict the analysis to a filling factor = 1/2, 
namely chose the chemical potential fi such that it results 
in half filling. 



D. Extension to a higher dimensional optical lattice 

We conclude this section by discussing the situation 
of a two dimensional optical lattice. The two fields are 
assumed to share the same wavelength and both inter- 
act with the same dipole transition of the atom. The 
Hamiltonian takes the form 



2D 



A 

2" 

2gx cos(i;) -|- 2gy cos(y) 



2gx cos{x) + 2gy cos(y) 
A 



(31) 

The lowest energy band in the symmetrical case with 
5a; = ffy = 0.1 and A — 0.25 is shown in Fig. [5l As 
for the one dimensional case, for q^, Qy — ±1/2 does the 
lowest energy band attain its maximum. Importantly, 
multiple minima of the dispersion is also found in the two 
dimensional case. It is easy to convince oneself that this 
holds also in three dimensional, and consequently that 
the topological FT is not limited to the one dimensional 
lattice. 




FIG. 7: First energy band Ei{qx,qy) of the two dimensional 
model ((31}, with equal coupling strengths = Qy = 0.1 and 
A = 0.25. 



in the collective atomic inversion. Moreover, it was found 
to be a non-zero temperature FT. 

To conclude, we showed that quantum PTs can oc- 
cur in multi-component ultracold atomic systems despite 
the lack of atom-atom scattering. It is believed that 
adding interaction between atoms can give rise 
to novel phenomena. In such a scenario, tv^ro 
relevant questions directly come to mind: How 
does the coupled internal structure of atoms af- 
fect the phase diagrams of non-coupled models, 
and how does strong interaction influence the 
topological PT demonstrated in this paper? The 
second question was briefly discussed at the end of 
Sec. lIII^ and both questions are currently under 
investigations for future publications. We further 
plan to analyze the present system, and especially the 
topological FT, in a cavity QED model where the field is 
treated quantum mechanically. 



IV. CONCLUSIONS 



In this paper have we studied an ideal gas of ultra- 
cold two-level atoms coupled to an optical lattice. It was 
shown that the coupled two-level character of the prob- 
lem give rise to novel phenomena, not present in the reg- 
ular dispersive case of large detuning where one atomic 
level has been adiabatically eliminated rendering an in- 
ternal structureless system. In particular, the Brillouin 
zone is twice as large in this system compared to a inter- 
nal structureless one and the Bloch bands possess mul- 
tiple minima. This comes about due to the competition 
between the various involved terms of the Hamiltonian; 
kinetic and internal energy. An outcome of the peculiar 
energy band spectrum is the presence of a topological, or 
Lifshitz, FT as the detuning is tuned across resonance. 
This transition is of first order nature and is manifested 



APPENDIX A: EFFECTIVE DISSIPATIONLESS 
TWO-LEVEL MODEL 

Decay of the excited atomic level will typically set lim- 
its on any experimental consideration. One way to min- 
imize such effects is to effectively couple two metastable 
atomic states using an additional laser. We consider 
therefore a three-level A-atom with metastable lower 
states |1) and |2) and excited state |3), with respective en- 
ergies Ei, i — 1, 2, 3. State 1 and 3 are coupled through 
our optical lattice, while 2 and 3 couple via an external 
laser whose field amplitude is assumed constant over the 
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atomic sample. The Hamiltonian becomes [28 



+2g cos{x + uJot) ((Ti3 + 0-31 ) 



(Al) 



where a, = |2)(2| - |1)(1| and = |1)(2| + |2)(1|. We 
have numerically confirmed that the band structure of 
(jA2p looks very similar to the ones presented in Fig. [1] 
Therefore, our results and conclusions of Sec. lIIII are also 
reproducible for a Hamiltonian such as (|A2[) . 



where f2 is the external laser coupling amplitude, lol 
and ojo the two field frequencies, and aij ~ |*)0|- The 
excited state can be adiabatically eliminated if the in- 
teraction is dispersive. Thus, we assume a highly de- 
tuned configuration, Ai — £'33 — E22 — lul ^ and 
A2 = £'33 — Ell — Wo ^ 2(7, and at the same time 
A3 = |Ai — A2I is of the same order as il and g. One 
then derives, after application of a rotating wave approx- 
imation, an effective two-level model for the metastable 
states, which is given by [11] 
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